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With the one-boson-exchange model, we study the interaction between the S-wave D M /D_* meson and S- 
wave B^'/Si*' meson considering the S-D mixing effect. Our calculation indicates that there may exist the 
6 t -like molecular states. We estimate their masses and list the possible decay modes of these B r -like molecular 
states, which may be useful to the future experimental search. 
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I. INTRODUCTION 



II. THE EFFECTIVE POTENTIAL OF B L -LIKE 
MOLECULES 
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Carrying out the study of the hadron configuration beyond 
the conventional qq meson and qqq baryon is an intriguing and 
important research topic. In the past decade, more and more 
charmonium-like or bottomonium-like states were observed 
in the e + e~ collision ILU-yll, B meson decays I0-01 and even yy 
fusion processes 081- UOll . which have stimulated the extensive 
discussion of exotic hadron configurations (for a review see 
Refs. EMI). 



In this work, we report on the investigation of 
hadronic molecules with both open charm and open bot- 
tom, where the interaction between the charmed meson 
(£)(*) = [D^,D^ + ,D ( * )+ ]) and bottom meson (S w = 
[fi w+ , BS*^, 5 ( / )0 ]) occurs via the one boson exchange (OBE). 
These new structures are labeled as the B f -like molecules be- 
cause such systems contain a charm (c) quark and an anti- 
bottom (b) quark. Because of the special hadron configu- 
ration, the prediction of the B e -like molecules with masses 
above 7 GeV can provide important information for further 
experimental search at facilities such as LHCb and the re- 
cently discussed Z° factory fll4fl . 



This paper is organized as follows. After the introduc- 
tion, we present the formulas of effective potential of B e -like 
molecules. In Sec. [TTT1 the numerical results are given. This 
work ends with the discussion and conclusion. 
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FIG. 1: (color online). The flavor wave functions of these hadronic 
molecular states, which consist of two isosinglets (a s i = (Z) W °B ( * )+ + 
DW + JSW°)/V2, a sl = £» ( ; )+ S ( ; )0 ), an isotriplet (a, = [D M °B M °, 
(D W °B W+ - £> W+ B W °)/ V5, D w+ B w+ ]) and two isodoublets (a<n = 

[fl W +B (,» D w +5 w +]; aj2 _ [ D (*)+b { ; )0 ,D mo B { ; )0 ]), where the index 



a is taken as X, Y, Z and Z corresponding to the T)S, !D*S* 
and DS" systems, respectively. 
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The B c -like molecules are categorized into four groups, i.e, 
£>S, D*!B\ D*S and £>£*. Each group contains nine states, 
which form an octet and a singlet. Their corresponding flavor 
wave functions are listed in Fig. Q] We adopt the approach de- 
veloped in Refs. il5l422ll to study the interaction of the B c -like 
molecules. In terms of the Breit approximation, the scattering 
amplitude iM{D ( * ) &* ) -> £) w £ w ) is related to the interac- 
tion potential in the momentum space by the relation 
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where M, and Mj are the masses of the initial and final states, 
respectively. The potential in the coordinate space 'V(r) reads 



as its Fourier transformation, 

VeQc) = J ^3 e* r T*(q) ^ 2 (<z\ m|) (1) 

where wig is the exchanged meson mass and the monopole 
form factor (FF) T{q 2 , m 2 E ) = (A 2 - m|)/(A 2 - q 2 ) is intro- 
duced to depict the structure effect of the vertex of the heavy 
mesons interacting with the light mesons. The parameter A, 
which is about one to several GeV, not only denotes the phe- 
nomenological cutoff, but also regulates the effective poten- 
tial. 

According to the heavy quark limit and chiral symme- 
try, the interactions of the light pesudoscalar, vector and 
scalar mesons with the S-wave heavy flavor mesons were con- 
structed as SQ 



(2) 
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(4) 



where the multiplet fields are expressed as H„ - -^-VP* a uY ' — 
Vajs], H^ = [P7y,-P a j5] l -4,H = yoWyo with v = (1,0), 



t)T 



(D^°,D^-,D ( s * y ) or (B<') + ,BW°,B ( ; )0 ), which satisfy the 



normalization relations <0|!P|g<7((r)> = (0\P\Qq(Q-)) = yjWp 
and {Q\P;\Qq(l-)) = (0\P;\Qq(l-)) = e, y(W>- The ax- 
ial current reads as A 11 = ±(£ f d M £ - %d^) = fd^f + ■■■ 

with € = exp(iP//„) and f n = 132 MeV. p^ = fe v V£,/ V5, 
Fpy(p) = d^pv-dvPn + lpn, p v ], F'^ip) = dfipv-dypn-lpn, p v ] 
and gy = trip/fw, with gy = 5.8. In the above expressions, P 
and V denote the three by three pseudoscalar and vector ma- 
trices, respectively, i.e., 
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The coupling constants involved in Eqs. (|2]i-(|4]i include 
0.59 extracted from the experimental width of D* + \5i 
[3 = 0.9 determined by the vector meson dominance mecha- 
nism, A = 0.56 GeV -1 obtained by comparing the form factor 
calculated by light cone sum rule with the one obtained by 
lattice QCD. In addition, the coupling constant related to the 
scalar meson cr, g s = g„/(2 Vo") with g n = 3.73 was given in 



Ref. 12811 . In the heavy quark limit, the interactions of the 
<£)(*)<£)(*) an( j <g(*)<g(*) v^ith light mesons are the same. 

With these Lagrangians listed in Eq. (f2]i-(|4]i, we can deduce 
the expressions of T^r). When obtaining the total effective 
potentials, we sandwich T^r) between the corresponding Be- 
like molecular states. Thus, the general expression of the total 
effective potential is expressed as 



^rl/W = ( a dJ]\ Yj V E (r)\af[J] 



E '= K, TJ,p,CO,, 
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where subscript a^ with £ = s\ , s2, t, dl , d2 and a - X,Y,Z,Z 
is introduced to distinguish the total effective potentials of the 
molecular systems defined in Fig. [TJ J denotes the total angu- 
lar momentum of system (7 = 0, J - 1,7 = 0,1,2 for the !D&, 
Tf'Bj'D'B* and TfB* systems respectively). The definitions 
of |a^[7]\ are 

|^[0]) = \DSCS )), 

|Z f [l]) = (|£>*S( 3 Si)),|£>*£( 3 D 1 ))) r , 

|z f [l]) = (|©£*( 3 Si)),|;D£*( 3 Di))) r ) 

|y f [0]) = (|D*S*( 3 5 )),|£)*S*( 5 /Jo))) r , 

|7 f [l]) = (\D*S*eS l )),\D*!B*eD l )),\D*!B*eD l ))) T , 
|7 f [2]) = (\D*£\ 5 S 2 )),\D*!B*( l D 2 )),\D*!B*eD 2 )), 

\D*S\ 5 D 2 ))f, 

with 



DS*/D*S(- + Lj) = J] cZ LmL K'-L,m L ), 
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D*S*( + Lj) = ^ CZ saM C\ZJ$C,L,m L ), 
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where C\ M , , Cl M , and C s '" s t , denote the Clebsch- 
Gordan coefficients. The polarization vector for the vector 
heavy flavor meson is written as e™ = +-We™ + ie 1 ") and 

e ( "' = el" . In the above expressions, 2s+1 Ly is applied to de- 
note the total spin S , angular momentum L, total angular mo- 
mentum / of the £)(*>g(*> systems, while L = S and L = D 
are introduced to distinguish S-wave and D-wave interactions. 
Because of the S-D mixing effect, the obtained total effective 
potentials of the TfB, D&* and TfB* molecular systems are 
in matrix form. The total effective potentials are composed of 
subpotentials as shown in Table. Q] 

The expressions of the subpotentials are 
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-g 2 Y(A, m IT , r), V? s = - -fS 2 g\ T(A, m v ,r), 



V»* = -^7(A ) m .,r)diag(l,l), 



■J) 'B 



V 



--^ v Y(A,m v ,r)diag(l,l), 



V» a = -g;nJ]Y(A,m a ,r), 



TABLE I: The relation of the total effective potential < V T f ta ,(r) and 
the subpotentials. Here, w is taken as 3 and -1 corresponding to the 
states marked by the subscripts si and t, respectively. Since the total 
effective potential of the DB* systems is the same as that of the T)"B 
systems, we only show the result for Tf"B. We use - to denote the 
case when the OBE potential does not exist since no suitable meson 
exchange is allowed for these systems. 
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with Y(A,m E ,r) = £- r (e- mEr - e - Ar ) - \^c" Ar , where we 
use superscripts £)Wg(*) to distinguish these subpotentials 
for the different systems while the introduced subscripts P 
and V denote the corresponding light pseudoscalar and vec- 
tor meson exchanges respectively. ni£ denotes the mass of 
exchange meson. Matrices ^l[J], S[J] and C[J] are listed be- 
low with ^[[0] = diag(l, 1), #[1] = diag( 1,1,1), Jl[2] = 
diag(l, 1,1,1), S[0] = diag(2,-l), fi[l] = diag(l, 1, -1), 

f o V5) 

S[2] = diag(-l, 1,1,-1), C[0] = , C[l] = 

I V2 2 J 
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In addition, the kinetic terms for the 2? ( *>£) ( *) systems are 
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and C[2] 
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where A = -j-J^r 2 ■£-, A 2 = A - 4- wii, m 2 and ;«3 are the 
reduced masses of the corresponding systems. 



III. NUMERICAL RESULT 

With the above preparation, in the following we illus- 
trate the numerical results for the B c -like molecular systems. 
In order to obtain the information of the bound-state solu- 
tions (binding energy and root-mean-square radius) of sys- 
tems listed in Fig. [1] we need to solve the coupled channel 
Schrodinger equation with the deduced effective potentials, 
which can answer whether these B t -like molecular states ex- 
ist or not. Here, we adopt FESSDE, a Fortran program for 
solving the coupled channel Schrodinger equation OU 13211 . to 
numerically obtain the binding energy and the corresponding 
root-mean-square radius. Additionally we also use a MAT- 
LAB package MATSCE JH] to do a cross-check. Usually 
the OBE potential is suitable to describe the interaction of a 
loosely bound state. Thus, we require the obtained binding 
energy in the range of ~ -20 MeV and the cutoff in the 
range of 1 ~ 5 GeV when presenting the result. 

In Tablellll we list the obtained typical values of the bound- 
state solution of these B f -like molecular systems, while the 
dependence of the results on A is given in Fig. [2] Among 
the 24 cases shown in Table. |TJ we find that there exist the 
bound-state solutions only for 17 states: 

1. D&: We find the bound-state solution only for the X s \ 
and X S 2 states. Both of these states are of the same 
quantum number, i.e., I(J P ) - 0(0 + ). The values of the 
cutoff A is close to 1 GeV for the X s i state. For the other 
isosinglet X&, the bound-state solution appears when 
taking A ~ 3.2 GeV. 

2. D*!B/D!B*: The bound-state solution exists only for the 
four isosinglets Z s \, Z s2 , Z s \ and Z s2 with 0(1 + ). Since 
the effective potentials of the TfB and D&* systems 
are the same, the dependence of the bound solutions on 
A for Z s i and Z i2 are almost similar to those of Z s \ and 
Zsi respectively (see Fig. |2), The small difference of 
the reduced masses also results in the difference of the 
typical values listed in Table. [II] when comparing the 
results of the states marked by the same subscript si or 
s2. 



TfB*\ For the D*S* systems, there are 15 states. 
Among them we find 1 1 states with bound-state solu- 
tions, which include the isosinglets Ffj , Y J ^°, Y^ 1 , 
Y J = l , Y J = 2 , Y J = 2 , isodoublets Y J =°, Y J =°, Y J = l , and 

si ' si ' s2 ' d\ ' dl ' d\ ' 

isodoublets y/ =0 ,y/ =1 . 



TABLE II: The typical values of the obtained bound-state solutions for the £)(*'£>(*' systems. Here, A, E, and r RMS are in units of GeV, MeV, 
and fm, respectively. 



System State A/E/r RM s State A/E/r RM s 



1.3/- 1.28/2.58 3.2/-2.03/1.99 

DS X sl 1.4/-6. 14/1.35 X s2 4.0/-10.87/0.94 

1.5/-13.91/0.98 4.8/-21.90/0.69 



1.3/-1.32/2.60 3.0/-0.76/3.12 

£>£* Z sl 1.4/-6.21/1.34 Z s2 3.5/-4.96/1.32 

1.5/- 14.02/0.97 4.0/- 11.06/0.93 



System State A/E/r R MS 



D*S 



1.3/-2.10/2.07 

Z s i 1.4/-7.92/1.19 

1.5/-16.69/0.89 

3.0/-1.83/2.04 

Z.,2 3.5/-7.55/1.08 

4.0/-15.08/0.80 



System 



D'S' 



State A/E/r RM& State A/£/r RMS State A/E/r RMS State A/E/r RMS State A/E/trms 



1.25/-2.70/2.00 1.9/- 1.66/2. 12 2.5/-2.23/1.82 3.3/-0.94/2.61 3.4/- 1.62/2.03 

1.30/-8.82/1.25 Y J s f 1.95/-6.65/1.13 Yf =0 2.7/-8.04/1.01 Y%° 3.4/-5.10/1.14 Y J d f 3.5/-6.53/1.03 
1.35/-19.47/0.94 2.00/-17.77/0.74 2.9/-18.64/0.69 3.5/-12.41/0.75 3.6/- 14.64/0.70 



1.25/-3.59/1.76 1.96/- 1.05/2.64 4.27/-3.20/1.53 4.98/- 1.08/2.46 

Yfi 1 1.30/-8.98/1.22 Y' s f 2.05/-6.92/1.13 Yf =l 4.54/- 10.85/0.87 Y J d f 4. .99 '/ '-1. 27/2.26 

1.35/-17.40/0.95 2.14/-18. 81/0.74 4.81/-24.74/0.60 5.00/- 1.49/2. 10 



0.96/- 1.28/2.59 

Y'f 1.05/-7.99/1.25 Y J f- 

1.14/-20.12/0.90 



2.0/-2.50/1.80 

2.1/-7.25/1.14 
2.2/- 14.63/0.87 



We use a hand-waving notation, i.e., five-star, four-star, 
three-star and two-star etc to mark the states in order to in- 
dicate that the bound-state solutions exist when the cutoff pa- 
rameter A corresponds to the different values: A < 1.5 GeV, 
1.5 < A < 2.5 GeV, 2.5 < A < 3.5 GeV, 3.5 < A < 5 GeV 
respectively. In this way we categorize these states according 
to the numerical results listed in Fig. |2]and Table HH(see Table 
[jlI]for more details). Usually the cutoff A is taken around 1 
GeV, which is a reasonable value, especially in the deuteron 
case. Thus, a five-star state implies that a loosely molecular 
state probably exists. The mass spectra of the BD, BD*, B*D, 
and B*D* molecular states with the {Qq}{Q (,) q} configuration 
were studied with the QCD sum rule approach [34], which 
correspond to the above six five-star B c -like molecular states 
obtained in this work. 



In the following, we will discuss the allowed decay modes 
of these predicted B c -like molecular states that may be help- 
ful to the future experimental search. All the five-star states 
Xji, Z sl , Z s i, Y j ^°, Y ] = l and Y ] = 2 are the isosinglet with 
subscript si. Their decay modes are listed in the 2nd-7th 
columns of Table. IIVI respectively. In addition, the decays 



of the four four-star states 



y y=o 

1 si ' 



,J=0 



a/=1 



and Y J ~ 2 are 



shown in the 8th-l 1th columns of Table. IIVI respectively. In 
Table. IIVI we also give the decay modes of the remaining 
five three-star states. In these decay channels, the B C {\P\) 



and B' C (\P\) mesons are the mixture of the l l P\ and l 3 Pi 
states [35S]: \B C {\P X )) = \B c (\ l P x ))cosG + |fi f (l 3 Pi)>sin<9, 
l^(li'i)) = -|B c (l 1 Pi)>sin0+ |B c (l 3 Pi))cos6». At present 
only BJ1 1 Sq) was observed with a mass m(B c (l l So)) = 6277 
MeV IB6I1 . We adopt the theoretical values from Ref. 0511 
when giving the decay channels of these B c -like molecular 



states, i.e., M 



M, 



B,.(2 3 S,) 



6882 MeV, M, 



6333 MeV, M Bc&So) = 6842 MeV, 
bai'Po) = 6699 MeV > m baipo = 



6743 MeV, M B[ilPl) = 6750 MeV and M Bc(P p 2) = 6761 MeV 
113511 - In obtaining these decay channels, we have only consid- 
ered the ground state of the light meson. 



IV. DISCUSSION AND CONCLUSION 

In short summary, we have studied the interaction between 
the S-wave D w /M* 5 meson and S-wave B^/B^ meson in the 
OBE model. With the obtained effective potentials, we predict 
the existence of many B c -like molecular states where we have 
already included the S-D mixing effect. Besides estimating 
their mass spectrum, we also list their decay modes. 

For comparison, we list the bound-state solution for Y^°, 
Y J ^ 1 and Y^ 2 when considering the one-pion-exchange 
(OPE) potential only in Table. [V] The one pion exchange 
force provides the main attraction in the formation of the Be- 
like molecular state, which is consistent with the observation 
in Ref. JH. 
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FIG. 2: (color online). The variation of the binding energy E and root-mean-square radius r RM s with A for the S (,) S ( *' system. Here E and 
fRins are in units of MeV and fm. The superscript 7 = 0, J = 1,7 = 2 denotes the total angular momentum J. 



TABLE III: Summary of the S c -like systems. 



State I(J P ) Remark 



DS 



State I(J P ) Remark State I{J P ) Remark 



X,i 0(0 + ) * * * * * Y J s f 0(0 + ) * * * • * Y J = l 0(1+) * * ** 

X s2 0(0 + ) * • * 



Y J =° 0(0 + ) • *** Y J = l 1(1+) ** 



Z„ 0(1 + )***** F/ =0 1(0 + ) **** Y J = l \(l + ) ** 

Z.,2 0(1 + ) *** Y J d f ±(0 + ) *** Y J = 2 0(2+) ***** 

Z sl 0(1 + )***** F^° |(0+)*** F£ 2 0(2+)**** 

Z s2 0(1 + ) *** F/f 1 0(1 + )***** 



For the other five-star states, we find that there does 
not exist the bound-state solution only considering the 
sigma meson exchange. Further, we notice that the p 
meson exchange plays a much more important role in 
the case of the other five-star states. With X s \ as an 
example, the typical values of its bound-state solutions 
with the p meson exchange alone are (A, E, r^Ms) - 
(1.5, -1.52, 2.43), (1.6, -4.60, 1.50), (1.7, -8.97, 1.13). The 
comparison of these results and those listed in Table HTlindeed 
indicates that the p meson exchange dominates the X sl , where 
E, A, and trms are in units of Ge V, GeV, and fm, respectively. 

With Ffj as an example, we also examined the sensitivity 
of the results to the coupling constant in the OPE case. When 
adoptin g g — 0.885, which is 1.5 times larger than g - 0.59 
in Ref. 13011 . we have to lower the A value in order to get the 
similar binding energy to that in the case of taking g - 0.59, 
i.e., 



that the narrow resonance Y(3872) lies very close to the DD* 
threshold, which was first observed in the Jj^n^if invariant 
mass spectrum of the B — > KJ/t//7T + n~ process |4J|. Similarly, 
the Zsi and Z s \ states can decay into the B C (1 3 S \)nn mode. 
The 7(3940) state was observed in B -> KJ/if/co |5j while 
F(4140) in B -» KJ/tfrcp @]. Similarly, the predicted (7/=°, 
Fff \ Y J = 2 ) or (Y£°, Y J =\ Y ] = 2 ) may be searched for in the 
^(l 3 ^ i)« or B c (l 3 S \)(f> modes, respectively. 

In the future, it will also be important to calculate the 
branching ratios of the different decay modes. Moreover, 
the investigation of the B f -like molecular states in other phe- 
nomenological models is also very interesting. Hopefully the 
investigations presented in this work will be useful to an ex- 
perimental search of them, which will be an interesting re- 
search topic. 
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Thus, the effect of varying the coupling constant on the bound- 
state solution can be compensated by changing the A value. 

Most of the predicted B e -like molecular states can decay 
into a B c meson plus light mesons. It is possible to find these 
states in the corresponding invariant mass spectrum. Recall 
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TABLE IV: The decay modes of the predicted B c like molecular states. Here, / denotes that the corresponding decay mode is allowed. 
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TABLE V: The bound-state solutions of Y J =°, Y J = [ and Y J = 2 only 

si ' si si J 

considering OPE potentials. Here, E, A, and r RMS are in units of 
GeV, GeV, and fm, respectively. 
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